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Statement of Significance 

QUESTION: Isotope concentrations in polar ice cores have long been used to 
estimate paleotemperatures. The validity of such a “paleothermometer” , how- 
ever, may be compromised by variations in the atmospheric circulation between 
very different climates. In particular, we need to know the following: to what 
extent do the relative contributions of different evaporative sources to polar pre- 
cipitation differ between climates? 

APPROACH: We used a general circulation model (GCM) fitted with tracer 
diagnostics to quantify, for both present-day and last glacial maximum (LGM) 
climates, the contributions of various oceanic moisture sources to Antarctic pre- 
cipitation. 

SIGNIFICANCE AND IMPLICATIONS OF FINDINGS: The higher latitu- 
dinal gradient of temperature that characterizes the LGM leads to a relatively 
higher contribution of tropical moisture to polar precipitation during that time. 
This would imply that LGM temperature reconstructions from ice cores may 
be biased slightly low. Further analysis is necessary, however, due to the uncer- 
tainty that exists in the LGM boundary conditions used to force the GCM. 


RELATIONSHIP TO MTPE SCIENCE PLAN 

Analysis of isotope/climate relationships is central to a reliable reconstruction of 
paleoclimatic fields, and establishing these fields is critical to an understanding 
of long-term natural climate variability and sensitivity. This work therefore 
addresses MTPE concerns regarding natural variability and change in long term 
climate. 
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Popular Summary: “The Origin of Antarctic Precipitation: A Modeling Ap- 

proach”, by G. Delaygue (CEA Saclay, France), V. Masson (CEA Saclay, France), 
J. Jouzel (CEA Saclay, France), R. Foster (Code 974, NASA/GSFC), and R. Healy 
(Woods Hole Oceanographic Institution). 

Isotope concentrations in polar ice cores have long been used to estimate paleotem- 
peratures. Underlying the use of this “isotope paleothermometer” is the assumption 
that the relationship between surface temperature and isotope concentration over 
time at a single geographical point is the same as that observed over space during the 
present-day climate. The validity of this assumption may in fact be compromised by 
several factors related to climate change. 

The specific factor studied in this paper involves the evaporative sources for polar 
precipitation. Climatic changes in the relative strengths of these sources would imply 
a need for a recalibration of the paleothermometer. To quantify such changes, we 
performed two GCM simulations, one of present-day climate and the other of the 
climate during the Last Glacial Maximum (LGM), roughly 18000 years ago. Evap- 
orative sources of Antarctic precipitation were established using special tracer diag- 
nostics. Results suggest that polar precipitation during the LGM does indeed consist 
of (relatively) more water from tropical oceans, a direct reflection of the LGM’s in- 
creased equator-to-pole temperature gradient and its increased sea ice extent, which 
reduces high latitude evaporation. This result implies that an uncalibrated ice core 
paleothermometer would produce LGM temperatures that are biased slightly low. 

Because LGM boundary conditions are still under debate, we performed a third 
GCM simulation using a modified set of LGM boundary conditions. Using this simu- 
lation gives some qualitatively similar results, though the tropical contribution is not 
quite as high. Uncertainties in the LGM boundary conditions does hamper success 
in calibrating the paleothermometer. 
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ABSTRACT (210 words) 


The contribution of different moisture sources to Antarctic precipitation for 
present-day and glacial conditions is estimated with the NASA/GISS Atmospheric 
General Circulation Model. Despite its low horizontal resolution (8°xl0°), this 
model simulates reasonably well the broad features of the observed present-day 
hydrological cycle. Simulated present-day Antarctic precipitation is dominated 
throughout the year by moisture from a subtropi:al/midlatitude band (30°S-60°S). 

The moisture supplied to a given coastal area of Antarctica originates mostly in 

* 

the adjacent oceanic basin; closer to the pole, other oceanic basins can also 
contribute significantly. 

Replacing the present-day sea surface temperatures (SSTs) and sea ice cover in 
the GCM with those from the CLIMAP oceanic reconstruction for the last glacial 
maximum (LGM) greatly increases the simulated latitudinal temperature gradient, 
with the consequence of slightly enhancing the contribution of low latitude 
moisture to Antarctic precipitation. It also changes the seasonality of the different 
contributions and thus their budget, particularly in coastal regions. Because the 
nature of LGM tropical SSTs is still under debate, we performed an additional 
LGM simulation in which the tropical SSTs are reduced relative to those of 
CLIMAP. The resulting decrease in the latitudinal gradient brings the relative 
contributions to Antarctic precipitation more in line with those of the present-day 
simulation. 



1. INTRODUCTION 


The Antarctic ice sheet, the largest surficial freshwater reservoir on the Earth 

(-30.10 6 km 3 ), is an important component of the climate system. Fluctuations of 
the hydrological cycle may induce changes in the ice sheet mass balance, with 
potentially important consequences for sea level change on short and long 
timescales (a sea level increase of about 70 m would result from a total melt ) 
(Warrick et al., 1996). Moreover, the Antarctic ice sheet constitutes a unique 
archive of past climate and environmental changes, with ice core records now 
extending through the last four climatic cycles (Petit et al., 1997) , roughly 
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400000 years. Figure .1 shows the various Antarctic sites where deep drillings 
have already been achieved (Byrd, Dome C, Dome B, Vostok and few coastal 
sites) , or is ongoing or planned (Dome Fuji, Siple Dome, Dome C and Dronning 
Maud Land; these last two drillings are undertaken under the framework of 
EPICA, the European Project for Ice Coring in Antarctica). 

The reconstruction of past Antarctic climate is largely based on the 
interpretation of the oxygen 18 (8 18 0) or deuterium (8D) profiles measured along 
these ice cores (Lorius et al., 1979; Jouzel et al., 1996). An empirical relationship 
between isotope concentration in precipitation and local temperature has been 
established from present-day measurements; this relationship has been used to 
interpret the isotope paleodata in the ice cores in terms of past temperatures. The 
implicit assumption in this approach is that the measured present-day spatial 
relationship holds at a single location through time. This assumption is now being 
challenged by alternative estimates of glacial to interglacial temperature changes 
inferred from borehole thermometry for Greenland (Cuffey et al., 1995; Johnsen 
et al., 1995) and to a lesser degree for Antarctica (Salamatin et al., 1998). 

Indeed, the assumption requires close scrutiny because the isotopic content of 
polar precipitation depends not only on the local temperature, but also on the 
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origin and transport history of the precipitating water (Koster et al., 1992). Thus, 
if the origin of precipitation changes in time, the temporal isotope-temperature 
relationship may differ from the present-day spatial relationship. Temporal 
changes in precipitation seasonality and vertical atmospheric structure may also 
lead to disagreement in the temporal and spatial relationships (see Jouzel et al., 
1997, for a recent review) , but Krinner et al. (1997a) have shown them to have 
little effect. 

A better understanding of the moisture sources for Antarctic precipitation - and 
how the relative contributions of these sources may have varied with time - would 
presumably improve the interpretation of the Antarctic ice core data. 

4 

Unfortunately, moisture contributions from different evaporative sources cannot 
be isolated through direct measurement. One common indirect approach is based 
on the relationship between evaporative source conditions (temperature, humidity, 
wind speed) and deuterium excess (d = 5D - 8 * 8 O) at the precipitation site. 
Petit et al. (1991) and Ciais et al. (1995) applied a Rayleigh-type isotopic model, 
which captures this relationship well, to Antarctic surface snow measurements and 
inferred a major subtropical evaporative source. They could not discount, 
however, an additional contribution from nearby coastal waters. More recently, 
using the same model, Delmotte (1997) suggested a seasonal variation in the 
relative strengths of the sources contributing to precipitation at the Law Dome 
site on the East Antarctic coast. Other studies concerning Antarctic precipitation 
origins, based on hydrological and isotopic approaches, give divergent 
conclusions (see Bromwich, 1988, for a review). Since then, an analysis of storm 
tracks across the southern hemisphere (Jones and Simmonds, 1993) has imposed 
new constraints on vapour origin, at least in a longitudinal sense. 

These indirect approaches generally assume an idealized moisture history 
(precipitation and mixing) along the trajectory from ocean to ice sheet. A more 
realistic representation of the moisture history can be obtained with an 
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atmospheric general circulation model (AGCM) fitted with special tracer 
diagnostics. With such a model, moisture evaporating from a specified source 
region can be "tagged", and this tracer vapour can be followed through the 
atmosphere until it precipitates. Examination of the tracer precipitation 
diagnostics at a given site establishes the evaporative sources for the precipitation 
there, as simulated by the AGCM (Joussaume et al., 1986; Koster et al., 1986). 

Using this approach, Koster et al. (1992) established the sources of 
Antarctic precipitation for present-day July conditions. They found that the 
average source temperature for July precipitation in Antarctica is roughly 12°C. 
Here, we extend this analysis to the full annual cycle and to both present-day and 
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ice age conditions in order to quantify the potential seasonal and climatic 
variations in the relative strengths of the moisture sources for Antarctic 
precipitation. 


2. TRACING WATER WITH AN AGCM 

2.1. The GISS atmospheric general circulation model: 

For the present study, we use the Model II version of the NASA Goddard 
Institute for Space Studies (GISS) atmospheric GCM (Hansen et al., 1983). The 
boundary conditions forcing the atmosphere dynamics include diumally-var- 
ying and seasonally-varying insolation, seasonally-varying sea surface 
temperature and sea ice distribution, and atmospheric composition (C0 2 content, 
etc.). Each grid box is divided into fractions of open sea, sea ice, land and 
continental ice, which are treated independently. Several studies have shown that 
this GCM, even when run at a coarse (10°x8°) resolution, successfully 
reproduces the main features of the earth's climate system, due in part to 
detailed treatments of physical processes. Our simulations use the 10°x8° 
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resolution, with nine o-levels in the vertical. We discuss in section 3 the 
simulated hydrological characteristics around Antarctica. 

2.2. Use of the tracer version of the GISS AGCM: 

Tracer diagnostics have been incorporated into the GISS AGCM to track the 
water evaporating from specified source regions as it is transported through the 
atmosphere (Koster et al., 1988). This is achieved by applying the same model 
physics to the tagged tracer water as to the model's standard water prognostic 
variables. This tracer model has also been used to transport water isotopes (Jouzel 

et al., 1987), using appropriate fractionation physics at each change of water 

€ 

phase. The resulting fields of isotope concentration in precipitation agree well 
with observed fields, lending some support to the simulated hydrological cycle. 

One potential problem in the application of the tracer model to the 
determination of Antarctic moisture sources involves the singularity at the south 
pole grid point in the GCM's cylindrical coorcinate system, since atmospheric 
transports across the pole is poorly represented. The very low temperature at the 
pole itself, however, reduces the air vapour content to half that at surrounding 
cells on average, so that any transport across the pole would, in any case, be 
small. 


In the present experiment, because the northern hemisphere's contribution 
to high southern latitude precipitation is fourd negligible (less than 1%, as 
determined from a similar simulation in which the northern hemisphere was 
marked as a single tracer), our evaporative moisture sources are defined only in 
the southern hemisphere. We first used basin geography to group the southern 
hemisphere's oceanic grid cells into the three basin-specific oceanic source 
regions shown in Figure 2a. All nonpolar cells ihat contain more than 70% land 
were combined into a fourth region, and grid cells within the Antarctic ice sheet 
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and the ocean immediately surrounding it (which is strongly influenced by sea 
ice) comprise two additional regions. The Antarctic ice sheet is isolated as a 
source region in order to assess the importance of polar moisture recycling. For 
some additional analysis, we regrouped the grid cells comprising the three ocean 
basin source regions into four source regions, defined by their annual SST. 
Figure 2b shows this set of source regions, labelled ‘equatorial’, ‘tropical’, 
‘subtropical’, and ‘midlatitude’. 


3. SIMULATION FOR MODERN CONDITIONS 

* 

In a 10-year present-day simulation, we forced the tracer GCM at the ocean 
boundary with the 1979-1988 SST dataset compiled by Reynolds (1988) , used in 
the Atmospheric Model Intercomparison Project (AMIP; Gates, 1992). The 
prescribed SSTs thus had realistic interannual variations. Some doubt being cast 
on the interannual sea ice variability, the 1979-1988 average was used instead for 
each month sea ice cover. The atmospheric C0 2 concentration was set to its early 
industrial level (315 ppmv). All other boundary conditions were as described in 
Hansen et al. (1983). 

As an indirect test of the model's ability to transport water tracers correctly, 
we now examine the resulting simulated water cycle around Antarctica in terms 
of (1) moisture transport and (2) snow accumulation. In Figure 3, the simulated 
zonally-averaged meridional transport of vapour in the southern hemisphere is 
compared with estimates from operational analyses of the ECMWF model 
(Bromwich et al., 1995). A previous attempt to estimate this flux from upper air 
station data (Peixbto and Oort, 1983) is also shown to emphasize the uncertainty 
due to difference in technique, though we note that the ECMWF analysis is 
considered far more reliable. The transport simulated by the tracer model is in 
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good agreement with the ECMWF estimates at nost latitudes, but it is too high 
over the Antarctic coast (between 65 °S and 75°vS). Indeed, at 70°S, the simulated 
poleward flux is 9.2 kg/m/s, twice as large as the estimate of 4.6 kg/m/s calculated 
by Giovinetto et al. (1997) from different numerical analyses. This excessive 
inland penetration of moisture is probably due to the coarse resolution of the 
tracer model, which allows an excessive numerical diffusion of moisture and does 
not simulate the barrier effect of the steep Antarctic coast (Genthon, 1994). 

Nevertheless, the model simulates reasonably well the large scale moisture 
transport from low latitudes to Antarctica. The longitudinal distribution of this 
transport is mainly controlled by the position of quasi-stationary cyclones along 
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the coast (Bromwich et al., 1995). Figure 4 shows the simulated annual 
meridional moisture flux across 70°S and the simulated annual sea level pressure 
around Antarctica, both averaged over the 10-year run. Although the absolute 
pressure values are too high by about 20 hPa, the pattern of pressure lows 
corresponds to the numerical analysis-based estimates of Bromwich et al. (1995). 
The stationary behaviour of these lows in the simulation is confirmed by a small 
associated variability; the standard deviation of sea level pressure is roughly 5 to 
10 hPa on both seasonal and interannual timescales. 

The meridional fluxes across 70°S are, as expected, well correlated with the 
local cyclonic pattern, i.e., they are systematically southward on the eastern side 
of a low center and northward on the wester i side. The strongest simulated 
moisture transport is associated with the Bellingshausen-Amundsen and Weddell 
Sea pressure lows, in agreement with other studies, although the Weddell Sea low 
and transport are too pronounced (e.g. Bromwich et al., 1995). It is not clear, 
however, whether these strong fluxes in nature result only from regional 
conditions (low pressure centres; latent heat flix from oceanic evaporation) or 
from cyclonic track feed (Jones and Simmonds, 1993). The model does shift 
westward by 30° the cyclonic centre that exists at about 120°W (Amundsen Sea), 
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and this could unrealistically decrease the contribution of East Pacific moisture 
into West Antarctica. Overall, though, the model simulates a correct longitudinal 
distribution of moisture inflow from the coast to the ice sheet. 

The simulated snow accumulation over Antarctica is 30 cm/year, which is about 
twice as large as estimates derived from operational numerical analyses (Budd et 
al., 1995; Bromwich et al., 1995) and AGCMs run at higher resolution (Genthon, 
1994). The excessive accumulation corresponds, in fact, to the excessive 
meridional vapour flux at 70°S, shown earlier in Figure 3. Our concerns about 
the excess annual accumulation are mitigated somewhat by the fact that the model 

reproduces a realistic seasonal cycle of precipitation (Figure 5). The model 
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simulates maxima of both precipitation and net accumulation in early fall (March- 
April), a precipitation minimum in spring (September) and an accumulation 
minimum in summer (December- January), which compares well with numerical 
analysis (Bromwich et al., 1995) and high resolution modelling studies (Tzeng et 
al., 1994; Krinner et al., 1997b). 

Figure 6 shows that in the areas surrounding the Vostok station in central 
Antarctica (represented by a 4 grid cell average) and the Law Dome station on the 
coast (2 grid cells), the magnitudes and seasonal cycles of the precipitation are 
reasonable when compared to the data of Dolgina and Petrova (1977) for Vostok 
and Casey/Wilkes (close to Law Dome). The simulated seasonal amplitudes are a 
little weak, but we note that the quoted observations rely, in part, on wind 
corrections to gauge measurements, which can sometimes increase the observed 
amplitudes dramatically. More recent observations for Law Dome (Vin Morgan, 
personal communication, 1997) suggest a significantly weaker seasonal cycle for 
the precipitation. 

The simulated accumulation shows an overly weak interannual variability. The 
standard deviation of annual accumulation is about 0.6 cm/year, which is much 
lower than the 30% interannual variation inferred from ECMWF analyses by 



Cullather et al. (1996). We do not expect that the inclusion of interannual 
variability in sea ice cover (neglected in the present simulation) would increase 
significantly the variability of accumulation, a similar low one was found in a 
simulation that includes these sea ice variations (Krinner et al., 1997b). 

In summary, the GCM simulation of the hydrologic cycle has some important 
biases that must be considered when interpreting the model results. Overall, 
though, model accuracy is deemed sufficient for the first-order analyses described 
below. 

3.1. Latitudinal contribution to precipitation: 

< 

The relative contributions of the 6 ‘latitudinal’ sources defined in Figure 2 to 
Antarctic precipitation are shown on Table 1. Also included are relevant physical 
characteristics of each evaporative source, namely the area, surface temperature, 
and local water budget. Notice that differences in the areal extent of the four 
warmest oceanic sources cannot explain more than 15% of their listed differences 
in evaporative contribution. Total evaporation from each source represents its 
potential role in providing vapour to the atmospheric water cycle. Evaporation 
minus precipitation indicates the net export of vapour to the atmosphere, but this 
cannot be interpreted in term of vapour origin because part of the local 
precipitation may originate elsewhere. 

The ‘equatorial’ zone provides less than 5% of the moisture precipitating in 
Antarctica. The ‘tropical’ zone, which is the strongest tracer source (evaporation 
is largest), contributes less than 15%. The small contributions of these two zones 
are consistent with the nature of the Hadley circulation, which keeps much of the 
tropical vapour away from southern latitudes. The small contribution of the 
‘equatorial’ region in particular can be explair ed by a high local turnover of 
moisture, as indicated by the small net moisture divergence (E-P) there. 



The ‘subtropical’ and ‘midlatitude’ zones (between roughly 30 and 63 °S) are 
the main contributors to Antarctic precipitation, providing up to two-thirds of the 
total. Notice that the ‘midlatitude’ zone, despite its large contribution, is 
effectively a vapour sink, with a negative net moisture divergence. This implies 
that much of the precipitation onto this oceanic region is simulated as derived 
from external (presumably equatorward) sources. In effect, much of the water 
travelling poleward from these external sources rains out, and some of it is 
replaced by the evaporated midlatitude water. 

Despite its low temperatures and correspondingly low evaporation rates, the 
polar oceanic zone around Antarctica supplies up to one fifth of Antarctica’s 
precipitation during austral summer and fall (Figure 9). The source contribution 
during austral winter is significantly lower. This seasonality corresponds well 
with that of sea ice area (e.g. Giovinetto et al., 1997, fig. 4) ; the evaporation 
being highest when sea ice area is smallest. Most of the local oceanic moisture is 
derived from the Weddell and Ross seas, where leads and polynyas 
(parameterized as fractional sea ice in the model) can allow an important air-sea 
exchange. 

Continental water recycling in Antarctica is not significant on the annual scale, 
as indicated by the small contribution listed in Table 1 for the "Antarctic" tracer. 
This small contribution stems from the extremely cold surface temperatures 
there, which lead to low sublimation rates. This contribution is significant (up to 
20% of the total) only during austral summer in central Antarctica, where the 
advection of moisture is small. 

One measure of the interannual variability of a relative contribution is the ratio 
of the standard deviation of the contribution to its ten-year average. In different 
Antarctic places, this ratio is generally less than 30%, which is similar to the ratio 
describing the interannual variability of measured Antarctic precipitation 
(Bromwich, 1988). Higher values of the ratio are usually found for the 



evaporative sources contributing the least; e.g., I he ratio for the Antarctic source 
is close to 70%. 

The contribution strengths discussed above are calculated for areally-averaged 
Antarctic precipitation. When tracer precipitation rates at individual grid cells are 
analyzed, we find in addition that the relative importance of the different 
evaporative sources varies geographically across Antarctica. 

Figure 7 shows, for austral winter and summer, the variation of source 
contribution with latitude over East Antarctica. Notice that the relative 
contributions of the midlatitude and high latitude oceanic sources decrease toward 
the pole, whereas those of the lower latitude sources slightly increase. This is 
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consistent with an idea earlier espoused by Petit et al. (1991), Koster et al. (1992), 
and Ciais et al. (1995) : moisture from different oceanic sources is generally 
transported at different atmospheric levels, and this allows the more remotely 
derived moisture to penetrate further poleward. 

Figure 7 also illustrates once again that the contributions from the different 
sources vary seasonally. This seasonality results from the interplay between 
cyclonic activity, which is intense in austral winter; sea ice growth, which 
modulates the vapour supply from the cold ocean, with a peak in fall; and 
continental water recycling, which is small but relevant in summer. 

Looking at the annual contributions to precipitation across Antarctica, we can 
compute for each grid cell an average source temperature, T-ave, by multiplying 
each source average temperature, Ti, by the soi ree contribution to the grid cell 
precipitation, Pi, and then dividing the sum of these products by the total tracer 
precipitation: T-ave=E(Ti.Pi)/E(Pi) . We consider only the oceanic sources for 
this calculation, since i) these sources supply mo* t of the precipitation water, ii) it 
allows a direct comparison with simple isotooic model, and iii) for isotope 
calculations, the relationship between source temperature and the isotopic 
character of the evaporate is stronger over the ocean. Figure 8 shows the spatial 



distribution of T-ave in Antarctica. These average source temperatures are 
slightly warmer than those calculated by Koster et al. (1992) with the same model 
(for July and with different SSTs), the averages of T-ave over Antarctica for July 
are 11.9°C (10-year std: 0.6 °C) in the present study and 11.6°C in the earlier 
one. The annual average is 11.3°C (with a 10-year standard deviation of 0.2 °C), 
reflecting the impact of seasonality on the source contributions: although sources 
are warmer during the summer, the coolest sources participate more so that the 
annual source temperature average is cooler than that for austral winter. Figure 8 
shows that the spatial distribution of T-ave is similar to that found by Koster et al. 
(1992). Slightly higher values of T-ave lie near the pole, reflecting a relatively 
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stronger contribution from remote sources there. 

These average source temperatures, by the way, are at odds with the results of 
Petit et al. (1991) and Ciais et al. (1995), who addressed the problem by analyzing 
snow isotope data in the context of a Rayleigh distillation model. Although they 
allowed for small contributions of moisture from higher latitudes and the coastal 
ocean, they inferred the dominant source of Antarctic precipitation to be the 
subtropical ocean, at 30°-40°S and 20°-40°S, respectively. Interestingly, using the 
same Rayleigh model, Delmotte (1997) inferred a significant contribution of high 
latitude sources to the precipitation at the coastal site of Law Dome. The Rayleigh 
distillation model, of course, is highly idealized and cannot account for the 
complexity inherent in atmospheric processes to the extent that a three- 
dimensional AGCM can. In particular Delmotte (1997) pointed out to the strong 
sensitivity of this simple model to isotopic initial conditions, poorly estimated in 
the previous studies (Jouzel and Koster, 1996). We thus consider our AGCM 
results to be more reliable. 



3.2. Basin contribution to precipitation: 

The relationship between each oceanic basin and climate variability is unique. 
The North Atlantic is often considered a source )f rapid climatic change through 
its control on the thermohaline circulation (e.g. Broecker, 1997). Data and model 
analysis (Blunier et al., 1998; Kim et al., 1998) , however, suggest that the ocean 
circulation in high southern latitudes may also play an important role. We might 
expect that climatic changes in zonal SSTs and meridional heat fluxes may differ 
significantly in the Atlantic, Pacific, and Indian Oceans and that such differences 
may be imprinted on the isotopic composition of the vapor evaporating from the 

basins. A major goal of EPICA is to examine the potential for such differences 
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through the analysis of long term isotopic profiles from two different sites: Dome 
Concordia, which is situated in the Indian Ocean sector, and Dronning Maud 
Land, in the Atlantic sector. 

Using the ocean basin tracers defined in Figure 2, we can examine with our 
present-day GCM simulation the contribution of each basin to Antarctic 
precipitation. Figure 9 shows these contributions, in percent, for the seasons 
exhibiting the greatest contrast. The GISS model clearly locates each basin domain 
of influence inside its longitudinal limits, the maximum influence being shifted 
eastward by 10 to 20° relative to the basin center, which is consistent with the 
cyclonic circulation that supplies moisture to Antarctica. Each basin contribution 
decreases rapidly poleward, the 50% isolines are generally not far from the coast. 
Moisture from the Indian Ocean, however, peretrates far inland during austral 
winter, and Pacific moisture largely dominates West Antarctica. 

Figure 10 shows, as a function of month the breakdown of the tracer 
contributions to precipitation in the regions surrounding Vostok on the Central 
Plateau and the EPICA drilling site in Dronning Maud Land (75-77°S; 0-10°E). 
Results shown are averages over several grid cells (see caption). The Vostok 
region is mainly under the influence of the Indian Ocean, which contributes 30% 



of the moisture during austral summer and 60% during austral winter. The high 
latitude ocean source is important during austral summer, supplying up to 25% of 
the moisture then. The model simulates similar contributions (not shown) for 
Dome Concordia, the other EPIC A drilling site (74.5°S; 123°E), but with a 
higher contribution from the Indian Ocean (75% in July). 

Precipitation in the Dronning Maud Land region is less dominated by a single 
ocean basin source. The closest ocean, the Atlantic, does supply most of the 
moisture (42% on average and up to 60% at the end of austral winter), but the 
Pacific Ocean supplies 20-30% of the moisture throughout the year, and the high 
latitude ocean source is strong in austral fall, when it provides 30% of the 
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moisture. Such a strong contribution from a high latitude source was inferred by 
Turner et al (1995) for a different coastal region (the Rothera station on the west 
coast of the Peninsula) from synoptic observations: about half of the depressions 
that led to precipitation at this station during a year of study had formed south of 
60°S. Seasonality in the high latitude source contribution is clearly seen at 
Dronning Maud Land, which lies near the coast; the maximum contribution is in 
austral fall, when the sea ice cover is smallest. On the Central Plateau (Vostok, 
Dome C), such seasonality is damped by the lower latitude contributions, which 
increase in austral fall due to a more intense cyclonic activity. 

Our simulation for present-day suggests that places characterised by low 
accumulation rates (where long ice records are possible) and dominated by 
moisture contributions from a single oceanic basin are restricted to the eastern 
end of the Central Plateau. The western end is at a crossroads of different 
influences. But simulations for the Last Glacial Maximum, in the following, show 
that the contribution of the most important source, at each place, increases with 
ice age conditions. 



4. SIMULATION FOR LAST GLACIAL MAXIMUM (LGM, 21 KA 
BP) CONDITIONS 


4.1. Description of experiments and simulated climates: 

As noted in the introduction, a better understanding of how evaporative source 
temperature changes with climate may improve our ability to interpret isotope 
paleodata. With this in mind, we performed two additional simulations with the 
tracer GCM, each one focusing on the climate of the LGM. 

The first simulation used the CLIMAP (1981) reconstruction of topography, sea 

surface temperatures and sea ice cover. The reconstructed Antarctic elevation was 

< 

generally increased by several hundred meters and up to a thousand meters in the 
west. An ablation of 200-500 meters took place along the eastern coast. C0 2 
concentration was lowered to 200 ppmv (Raynaud et al., 1993) and orbital 
parameters were set to values corresponding to 21 ka (thousands of years) ago 
(Berger, 1978). Aside from this last change, by the way, the climatic aspects of 
this simulation are identical to those in the ice age simulation described by Rind 
(1987), which used insolation parameters corresponding to 18 ka years ago. 

The second LGM simulation was performed to address uncertainty in the 
reconstructed intertropical SSTs provided by CLIMAP — some (e.g. Rind and 
Peteet, 1985; Pinot et al., submitted) believe that these SSTs are significantly too 
high. The nature of the intertropical SSTs is important in our analysis, of course, 
given the sensitivity of the global hydrologic cycle to the spatial distribution of 
SSTs and particularly to their meridional gradient. In the second LGM 
simulation, we cooled intertropical SSTs by 2 to 3° relative to modem conditions: 
(a) in the Atlantic Ocean we applied a 3° cooling between about 40°S to 25°N; (b) 
in the Pacific Ocean we applied a 2° cooling from about 48°S to 38°S; and (c) in 
the Indian Ocean, we applied a 2° cooling in each cell north of 26°S. This 
somehow arbitrary reconstruction is based on independent proxies, such as 



alkenones and corals (Edouard Bard, pers. com., 1997), and it is roughly 
consistent with output from coupled atmosphere-ocean simulations of LGM 
climate (e.g. Ganopolski et al., 1998). This imposed cooling lies between the 
severe cooling seen on the continent and the sea surface stability that can be 
inferred from faunal assemblages. 

The prescribed conditions in both LGM simulations have no interannual 
variability, so both simulations were integrated for only 5 years, after a one 
month spin-up period. The analysis below focuses on the monthly averages 
derived from these 5 years. The interannual variability of the simulated fields is 

found to be comparable to that in the present-day run, which uses interannually 

* 

varying SSTs. Hence, we infer that glacial conditions increase the intrinsic 
atmospheric variability, partly through the transient eddy circulation (Rind, 
1987). 

The meridional vapour flux simulated with standard CLIMAP conditions is 
very similar to that simulated in the present-day run (Figure 11), except for a 10° 
southward shift of the transport maximum consistent with the largest latitudinal 
gradient of temperature (Rind, 1987). Compared to the CLIMAP simulation, the 
cooler tropics imposed in the second LGM simulation lead to lower high latitude 
temperatures (by 1-2°) , but the net effect of this cooling in low and high latitudes 
is a weaker meridional temperature gradient. The southward flow is reduced in 
this "cooler tropics" simulation (by 15%) compared to present-day, due to the 
decrease of both global evaporation and meridional transport. 

In both LGM simulations, the net accumulation over Antarctica is reduced by 
about 20% — the 50% reduction in annual sublimation induced by lower surface 
temperatures is not enough to overcome a concurrent reduction in precipitation (- 
22% for the CLIMAP simulation, -28% for the cool tropics one). Due to a drop 
in sea level, however, Antarctica’s area as reconstructed for the LGM is larger by 
18%. Averaged over this larger area, reductions in both precipitation and 



evaporation are smaller (-6% and -22% respeclively for CLIMAP conditions, - 
20% and -34% with cooler tropics). The seasonal variation of precipitation is 
shifted into phase with the solar cycle, with extrema in summer and winter 
concurrently with sublimation. As a result, the seasonal cycle of accumulation is 
quite muted. This seasonal shift of precipitation, which is also seen in the tracer 
contributions, may reflect our seasonal interpolation of the CLIMAP SSTs and sea 
ice cover, which assumes extrema in February and August and a sinusoidal cycle 
in between. 

In the CLIMAP simulation, the simulated cooling over Central Antarctica 
(-3.5° annually) appears to cause the simulated decrease in precipitation there. We 

4 

conclude this from an analysis in which we scaled the LGM precipitation rates by 
a factor equal to the ratio of the present-day saturated vapour pressure (as infer- 
red from surface air temperature and pressure) to that of the LGM. These scaled 
precipitation rates are quite similar to the simulated present-day rates in Central 
Antarctica, with a regression slope of 0.95 (r 2 =0.85, n=575). Such a correlation 
between surface temperature and precipitation in Central Antarctica was already 
inferred from estimates of past accumulation ard temperature in deep ice cores 
(e.g. Yiou et al., 1985). This correlation, however, does not hold over all of 
Antarctica; simulated LGM precipitation in West Antarctica is higher than 
expected from the cooling alone. This might l)e explained by higher cyclonic 
activity in the west and its provision of greater amounts of moisture there. 
Indeed, in an LGM simulation with a stretched-grid AGCM, Krinner et al. 
(1997b) simulated a smaller decrease of this eye Ionic activity along the coast of 
West Antarctica than anywhere else on the contir ent. 

4.2. Contributions from tracer source regions: 

For consistency, we used the geographical definition of the moisture sources 
shown in Figure 2 in the two LGM simulations, despite the fact that the SSTs 



within these regions are less homogeneous during the LGM. Results are shown in 
Table 1 for the whole Antarctic continent. The relative contribution of each 
source is similar to that found for present-day conditions. The similarity in the 
relative source contributions between climates also holds at each grid cell in 
Antarctica. Small differences between the runs are generally comparable in 
magnitude to the interannual variability and thus must be considered with care. 
Still, they form a consistent picture when considered along with other climatic 
differences. 

In the CLIMAP simulation, the moisture contribution from low latitudes is 
generally increased, particularly from the Pacific in West Antarctica. This is 

4 

perhaps related to the CLIMAP SST reconstruction in the tropics, which: (1) 
increases tropical evaporation, primarily through a warming of the Pacific, and 
(2) increases wind strength in the Atlantic through a larger meridional gradient of 
SST (Miller and Russell, 1990). Outside of the tropics, the simulated increase in 
the atmospheric transport of energy is quite limited (Rind, 1987). 

In the CLIMAP simulation with cooler tropics, the meridional temperature 
gradient is more similar to that in the present-day simulation, and the relative 
contributions from the source regions are in better agreement with those for the 
present-day climate. Still, though, the LGM gradient is a little higher than the 
present-day gradient due to a large cooling at high latitudes. The result is a 
slightly higher relative contribution of subtropical moisture to Antarctic 
precipitation. Note that in this LGM simulation, the evaporation itself is reduced 
at all latitudes. 

Table 1 details the source contributions to total Antarctic precipitation. At the 
regional scale, when going from present-day to LGM conditions, the changes in 
the relative contributions are much larger than implied in the table. For example, 
in the Central Plateau, in the regions surrounding Vostok and Dome Concordia, 
the contribution of the Pacific source doubles under CLIMAP conditions at the 
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expense of the high latitude ocean and Indian Ocean sources. This is probably due 
to i) the increase in sea ice cover at high latitudes, and ii) the warming of the 
tropical Pacific associated with a stronger meridional temperature gradient. In 
contrast, though it also shows a weakening of the high latitude ocean contribution, 
the ‘cool tropics’ LGM simulation features a dramatic increase in the Indian 
Ocean contribution (up to 60% on annual average, and never less than 40%). In 
Dronning Maud Land, with CLIMAP SSTs, the decrease of the high latitude 
moisture supply slightly enhances the Pacific and Atlantic proportions, by about 
5% annually. The ‘cooler tropics’ simulation decreases the Pacific contribution 

and increases that of the Atlantic to about 50% of the annual precipitation. 

* 

The annual weighted-average oceanic source temperatures, as defined in section 
3.1, are shown for both LGM simulations in Figure 8. Differences between the 
LGM and present-day T-ave values are also shown. As expected, the CLIMAP 
reconstruction increases the average source temperatures. The average of T-ave 
across the continent (weighted by precipitation amount) is 13.4°C, which is 2.1 ° 
higher than that for the present-day simulation. On the other hand, the ‘cooler 
tropics’ LGM simulation shows slightly lower values of T-ave, with a continental 
average of 11.1 °C. These differences reflect the complex interplay between two 
separate factors controlling the changes in T-ave: the general decrease in SSTs, 
which tends to decrease T-ave, and the greater importance of more remote 
sources, which tends to increase it. For both LGM simulations, the highest 
increase in T-ave is seen at the strongest stationary low location, in Ross Sea and 
off the east coast (see Figure 4), secondly in Weddell Sea. In essence, this increase 
results from the replacement of cold high latitude moisture by warm subtropical 
moisture. These regions of highest increase of T- ive also show the warmest T-ave 
for LGM. It suggests that these cyclonic regions ire especially able during glacial 
climates to transfer moisture from low to high latitudes through the Circumpolar 
Trough. 



5. SUMMARY AND DISCUSSION 


We have used an atmospheric general circulation model (AGCM) fitted 
with water tracer diagnostics to quantify the sources of Antarctic precipitation 
for both present-day and ice-age climates. Before summarizing the results, two 
important caveats must be underlined. First, this AGCM has a low horizontal 
resolution and was not built explicitly for polar studies. Our analysis could 
perhaps be improved through the use of a stretched-grid model (a 'zoom' version 
of the AGCM) and by taking better account of the complex interactions known to 
exist between the atmosphere, sea ice, and ocean. Second, the realism of our 
tracer results can only be estimated indirectly through a comparison of the 
overall simulated hydrology with observed hydrological variables. Fortunately, 
the basic character of the simulated water cycle in the southern hemisphere and 
Antarctica is in many ways similar to that ascertained from data and numerical 
analyses (e.g. Jones and Simmonds, 1993; Bromwich et al., 1995), and it is 
deemed adequate for the first-order analyses performed here. 

For the present-day climate, most of the moisture that precipitates in 
Antarctica originates in the ‘subtropical’ and ‘midlatitude’ zones defined in Figure 
2. The relatively small contribution of equatorial and tropical moisture to this 
precipitation mirrors similar results found for Greenland with this model (Koster 
et al., 1986; Charles et al., 1994). When the Last Glacial Maximum (LGM) 
climate is simulated with two different SST reconstructions, we find that the 
increased latitudinal temperature gradient and the increased sea ice extent result 
in a slightly higher contribution of low latitude moisture to Antarctic 
precipitation. This is particularly true in the CLIMAP simulation, which includes 
a warm intertropical zone (especially in the Pacific) that provides significant 
moisture. Cooling this intertropical zone in the second LGM simulation 
intensifies the cooling in high latitudes (by 1 to 2°) and decreases slightly the 
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meridional gradient, so the increase in the low latitude contribution is not quite as 
strong. 

The breakdown of source contributions in Table 1 is reflected in the 
average oceanic source temperature (T-ave) for annual Antarctic precipitation, 
computed to be 11.3°C for the present-day climate. An average source 
temperature can also be computed at each Antarctic grid cell. Figure 8 shows that 
grid cells closer to the coast have a slightly higher tendency to derive moisture 
from colder (less remote) ocean surfaces. The figure also shows that the average 
source temperature can vary significantly with climate. 

Our AGCM results show that in present-day and even more in LGM climates, a 

< 

given region within Antarctica derives much of its moisture from the proximal 
ocean basin, especially if the region lies near the coast. In the eastern part of the 
Plateau, the Dome Concordia site is dominated by Indian Ocean moisture, with a 
weak seasonality in the Indian Ocean's contribution stemming from coastal and 
local (ice sheet) contributions. In the western pait of the Plateau, Dronning Maud 
Land receives moisture mainly from the Atlantic and the Pacific Oceans, with a 
contribution from high latitude ocean that varies seasonally with the presence of 
ice cover. Note that potential model errors increase toward the pole, so that 
results there must be considered with more cauti on. 

The next step in our research is to examine how the seasonal and climatic 
variations in the contributions of the moisture sources relate to isotopic variations 
in ice core measurements. The AGCM simulations discussed above suggest some 
important clues and issues. In particular, the different values of T-ave computed 
for the two LGM simulations (Figure 8) suggest that the nature of intertropical 
SSTs during the LGM must be firmly established if the interpretation of isotope 
paleodata is to account properly for changes in source temperature. 
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Figure captions 

Figure 1 : The Antarctic ice core drilling sites and regions discussed in the text. 

Figure 2 : Definition of the moisture sources for the simulations. Between the 
Equator and 63°S, the definition is based on either basin geography or annual sea 
surface temperature. The ‘high latitude ocean’ source is characterised by a 
seasonal sea ice cover. 

Figure 3 : Annual zonally-averaged meridional transport of moisture by the 
atmosphere: comparison of estimates based on ECMWF numerical analyses by 

t 

Bromwich et al (1995), upper air station data by Peixoto and Oort (1983), and 
present-day simulation results with the GISS 8°X10° AGCM. 

Figure 4 : 1979-1988 annual averages simulated by the GISS model of (top) 
meridional northward moisture transport across 70°S latitude (in kg/m/s); and 
(bottom) sea level pressure (isobars every 2 hPa from 998 to 1010). Note the 
pattern of quasi-stationary lows around Antarctica. 

Figure 5 : 1979-1988 averages of precipitation, sublimation and net accumulation 
(in mm/month) simulated over Antarctica with the GISS AGCM. 

Figure 6 : 1979-1988 average of precipitation (in mm/month) for central and 

coastal Antarctic stations: simulated values for Vostok and Law Dome (weighted 
averages of 4 and 2 grid points, respectively) versus estimates of Dolgina and 
Petrova (1977) for Vostok and Casey/Wilkes. 

Table 1 : Average contributions of the differer.t zones (Figure 2) to Antarctic 
precipitation, as simulated by the GISS model fo * modem and ice age conditions. 
Annual average and seasonal extrema are given for each zone, as well as the zone 
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area (in 10 12 m 2 ), the spatial range of annual surface temperature (in °C), and the 

hydrological balance (fluxes in 10 15 kg/yr). The ice age simulations consist of a 
full CLIMAP simulation and a simulation in which cooler tropics are imposed. 

Figure 7 : Latitudinal variation of moisture source contributions to East 

Antarctica precipitation (1979-1988 average), as simulated by the GISS model for 
austral winter (JJA) and summer (DJF). Note the significant seasonal variation in 
the high latitude and Antarctic contributions. 

Figure 8 : Average oceanic source temperature, T-ave, in °C, for annual 

precipitation under present-day and LGM conditions. Differences in T-ave 
between LGM and present-day conditions are also shown. 

Figure 9 : 1979-1988 average of the contributions of different ocean basins (see 
Figure 2) to the precipitation simulated by the GISS model, shown as a percent of 
total precipitation. The two seasons with the most extreme differences are shown 
for each ocean basin. Isolines are at every 10%, from 10 to 60%. 

Figure 10 : Annual source contributions to two areas: the Vostok region (on the 
Central Plateau, with a simulated annual temperature of -49°C) and the potential 
drilling site of EPIC A in Dronning Maud Land (75-77°S, 0-1 0°E ; simulated 
annual temperature: -40°C). Results shown are weighted averages over 4 and 2 
grid points, respectively. Each error bar represents one standard deviation, as 
determined from the 10 simulation years. Contributions from the source missing 
in the figure (nonpolar continents) would bring the total to 100%. 

Figure 11 : Annual zonally-averaged northward flux of vapour (in kg/m/s) 
simulated by the GISS model for three different climates: present-day, 



representing the 1979-1988 period, the LGM using the CLIMAP reconstruction, 
and the LGM using a CLIMAP with ‘cooler tropics’ reconstruction (see text for 
details). 
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